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o : 

04 ' We show that ten(one) years operation of the ultimate DECIGO (DECihertz Interferom- 

^ I eter Gravitational wave Observatory) can determine the cosmic equation of the state with 

Q \ such accuracy that 0.06%(3%), 0.08%(4%) and 0.06%(3%) for Qm, Om and w, respectively. 

^\ • In more realistic case of practical DECIGO or BBO (Big Bang Observer), w will be de- 

'"^ ' termined within ~ 10% by ten years observation assuming the flat universe model. Hence, 

I the DECIGO or BBO will give an independent determination of the cosmic equation of the 

, state. 

(N ■ 
> ■ 

' §1. Introduction 

^ : 

\ The laser interferometers begin to search for the astrophysical gravitational wave 

' sources in broad frequency band. For the higher frequency band of 10 — 10^ Hz, the 

■ ground-based interferometers such as LIGO, TAMA, VIRGO, GEO are currently 
in operation, and the next generation of interferometers such as advanced LIGO 
and LCGT^) are planed. For the lower frequency band of 10""^ — 10~^ Hz, the 
spaced-based interferometer LISA (the Laser Interferometer Space Antenna) will 

O \ be launched in ~ 2012. To fill the frequency gap between the ground-based detectors 

and LISA of 10~^ — 10 Hz, the deci hertz laser interferometer in space is planed to be 
03 ■ constructed in ~ 2020. This detector is called DECIGO (DECihertz Interferometer 

Gravitational wave Observatory) in Japanese group^^'^^ or BBO (Big Bang Observer) 
in the NASA SEU 2003 Roadmap "Beyond Einstein" .5) 

The major scientific objectives of the decihertz antenna are: 
\ • The primordial gravitational wave background will be observed, since there are 

no or little confusion noise from galactic white dwarf binaries around 0.1 Hz.'^^'^^ 

• ~ 10^ chirp signals of coalescing binary neutron stars and stellar mass black 
holes per year will be detected. By analysing the signals from these binaries at 
cosmological distances, it may be possible to determine the expansion rate of 
the Universe and the equation of state for the dark energy.'^) 

• The merger of intermediate mass black holes (10^ — IQ^Mq) will be observed 
and hence it is helpful to understand the formation and growth history of super 
massive black holes. 

In this paper, we investigate the second objective, i.e. the determination of the 
expansion rate of the Universe by the decihertz gravitational wave detector. 

Recently, from the observation of distant supernovae the expansion of the Uni- 
verse appears to be accelerating at present, which provides the evidence for a non 
zero cosmological constant or more generally the dark energy with the negative 
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pressure (see Ref.8) and references therein). In near future, SNAP (the Super- 
nova/Acceleration Probe) will observe ~ 2000 supernovae per year for z < 1.7 and 
will determine the energy density of the dark energy and its equation of the state 
with high accuracy (within a few % error). ^) But the measurement of supernovae 
might have possible systematic uncertainties such as the dust extinction, the gravi- 
tational lensing and the metallicity dependence. Hence it is desirable to determine 
the nature of the dark energy by an independent method, since the nature of the 
dark energy is by far the most important to clarify its origin theoretically. 

Seto, Kawamura and Nakamura (2001) suggested that the acceleration rate of 
the Universe will be directly determined by a 10 yr observation of the gravitational 
wave from a neutron star binary at z = 1. This is because the phase of the signal from 
the distant source is modulated due to the cosmic acceleration and this modulation 
can be measured by using the matched filtering techniques. The host galaxy of the 
source can be identified by the follow-up observations using electromagnetic wave 
antennae since the angular resolution of the ultimate DECIGO is a few arcsecond 
even for the source with z ~ 1, and hence the source redshift is also determined. ^^*-' 
Thus the expansion rate of the Universe (i.e. Hubble parameter H{z)) is directly 
obtained at the source redshift z. 

Following Seto et al.(2001), we evaluate the accuracy of determining the cosmic 
acceleration rate by the decihertz gravitational wave antenna. We newly include 
the effects of (1) the number of NS/NS merger events and (2) the various masses 
and redshifts of the binaries. From the 1st effect, since ~ 10^ signals of NS/NS 
binaries will be measured per year, the accuracy of the determination of the cosmic 
acceleration is l/VTo^ ~ 10~^ — 10~^ times better than that for a single event. From 
the 2nd effect, we investigate the dependence of the accuracy on the binary masses 
(0.1 - W^Mq) and its redshifts (z = - 2). 

We adopt {{2rn,^A) = (0.3,0.7) cosmology with the Hubble parameter Hq = 
70kms~^Mpc~^ and use the units of c = G = 1. 

§2. Basics 

2.1. Sensitivity of DECIGO/BBO 

DECIGO/BBO would consist of three spacecrafts orbiting the sun, in a trian- 
gular configuration with separations of 5 x 10^ km, which is 100 times smaller than 
that of LISA.^) The strain sensitivity of the practical DECIGO or BBO is about 
1000 times better than that of LISA (i.e., 10^^^Hz~^/^), and acceleration noise is 
100 times smaller than that of LISA (i.e. 3 x 10~^'^m s~^ Hz"^''^). Using the above 
assumptions and the sensitivity curve of LISA in Ref.ll), we plot the sensitivity of 
practical DECIGO and BBO in Figdi^) We also show that of LIGO II (or LCGT). 
The sensitivity of the ultimate DECIGO is 1000 times better than that of the prac- 
tical DECIGO (i.e. h ~ 10~^^ at / = 0.1 Hz), which is determined by the quantum 



*' In Ref. 4), the source is at 300Mpc but the sensitivity is ~ 1000 times worse than the ultimate 
DECIGO. It is possible to predict the source position within an arcminute in Ref.*-* so that for the 
source at ~ lOGpc {z ~ 1) the ultimate DECIGO can determine the position within a few arcsecond. 
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limit sensitivity for a 100 kg mass mirror. ^-^ In this paper, we use the sensitivity 
curve of ultimate DECIGO, to show what we can do ultimately. 

2.2. Gravitational Waveform of In-spiralling Binary in Accelerating Universe 

We consider a binary of mass Mi^2 
at the redshift z as the gravitational 
wave source. As the binary emits 
the gravitational radiation, its orbital 
separation decreases to zero. From 
the Newtonian quadrupole approxima- 
tion,^^) the time left till the coales- 
cence is given by. At = t^ — t = 
5{8TTfy^/^Mz^^^, where t^ is a coa- 
lescence time, / is the frequency of 
the gravitational waves, and Mz = 
(MiM2)^/^(Mi-hM2)-^/^(l-Fz) is a red- 
shifted chirp mass. 

The time to the coalescence At at 
present is related to that (Atz) at the 
redshift z as,^)'^^) 
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Fig. 1. Sensitivity of DECIGO/BBO. 



AT = At- X{z)At'^ 



(2-1) 



where AT = Atz{l + z) and the second term is the correction term due to the cosmic 
acceleration. The acceleration parameter X{z) is defined by,^) 



X(z) 



Ha 



Hjz) 
l + z 



(2-2) 



where H{z) is the Hubble parameter at z. X is rewritten as X = (a(0) — a{z))/2 
where a is a scale factor and ' = d/dt. Thus if X{z) > [X[z) < 0], d(0) > d{z) 
[a(0) < d{z)] and the universe is accelerating (decelerating) (see Ref.3)). For the 
small redshift z <S 1, X{z) is simply reduced to X{z)/Hq = —{l/2)qQZ where go is 
the deceleration parameter. The phase of the gravitational wave signal at present is 
shifted by -Itt f X (z) At'^ from Ea.lQTl). 

In the frequency domain, the gravitational waveform including the effects of the 
cosmic acceleration is written as,^) 



kccif) = Hf)e'"'-^^\ 
where the phase correction is d^a_cc = —2'KfX{z)At'^ or, 

257r-i3/3 



32768 



-X(z)/"13/3a^;10/3_ 



(2-3) 



(24) 



h{f) in Ea. (|2-3() is the in-spiral waveform and we use the restricted first post- 
Newtonian approximation. Then, we have h{f) = Af~'^^^e'^'^^-f\ where A is the 
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amplitude and ^{f) is the phase. ^'^^ They depend on five parameters: the redshifted 
chirp mass Mz, the reduced mass /i^ = MiM2(l + z)/{Mi + M2), the coalescence 
time tc, the phase (pc, and the luminosity distance to the source -D5. The amplitude 
is given by ^ = K (5/96)^^'^ Ai^^^ / {tt'^/^ Ds) , where K is the constant determined by 
the inclination of the source, the relative orientation of the source, and the detector. 
Since the average value of K is about unity, we assume K = 1 for the following 
calculation. The phase ^{f) is a rather complicated function of M.z, ^J'z, (pc, and t, 
which is given in Eq.(3.13) of Cutler & Flanagan (1994). 



2.3. Parameter Estimation 

The signal haccif) in Ea. H2-3|) is characterized by six parameters 7^ = {Aiz, fJ'z, 
(j)c, tc, Ds, X). In the matched filter analysis with the template, these parameters can 
be determined. We compute the errors in the estimation of these parameters using 
the Fisher matrix formalism^^^' (see also Refs.l7)-21)). The variance-covariance 
matrix of the parameter estimation error Aji is given by the inverse of the Fisher 
information matrix Fij as {A^jiAjj) = {^r~^^.y The Fisher matrix becomes 

where Sn{f) is the noise spectrum. We regard Sn{f) as the instrumental noise of 
ultimate DECIGO in Fig^ neglecting the binary confusion noise since there is no 
or little confusion noise for / ^ 0.1 Hz.'^^ The signal to noise ratio {S/N) is given 
by 

{S/N)^ = Aj \hacc{f)^. (2-6) 

We integrate the gravitational waveform in Ea. (|2-5() and Ea. H2-6|) from 1, 3, and 10 
yr before the final merging to the cut-off frequency fcut when the binary separation 
becomes r = 6(Mi + M2). 

The estimation error A^y simply scales as {S/N)^^ from Ea. H2-5|) and Ea. ()2-6|) . 
Hence the errors for the practical DECIGO or BBO is simply 1000 times larger than 
that for the ultimate DECIGO (see FigHJ. 

2.4. Event Rate of NS/NS Coalescence 

The estimation error A'y discussed in the previous section is for a single event. 
But for a number of events the estimation error decreases in proportion to (the 
number of events)"^/^. For this purpose we first estimate the event rate of NS/NS 
coalescence. 

The NS/NS merger rate between the redshift (z, z + dz) is written as, 

dN , r{zf 

where r{z) is the comoving distance and ng is the comoving number density of 
the NS/NS merger rate per year. Here, we assume that the comoving merger rate 
no is constant. We choose tiq so that 1 event per year could be observed within 
60 Mpc from recent studies, ^^^'^^^ and hence we set no = 1 x 10~^Mpc~^yr^-'^. 
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In FigEJ the merger rate per redshift bin 
Az = 0.1, AN = {dN/dz)Az, is shown. 
Since the ultimate DECIGO could ob- 
serve the NS/NS binaries at cosmologi- 
cal distances {z ~ 1) with high S/N ~ 
10^, the NS/NS merger will be observed 
10^~^ events per year from Fig|21 Thus, 
the estimation error A^y is 1/^/AN (~ 
10"^ or 10"'^) times better than that for 
a single event. 

§3. Results 



^ We show the estimation error of 

x^- „ rr., ^.r, ,^rr, r thc cosmlc accelcratlou parameter Xiz) 

bis. 2. the JNb/Nb merger rate per year tor „ 

zi2 = 0.1bin. ^ ^ ^ in Eq.(|13 at z = - 2. We con- 

sider the mergers of NS/NS binaries 
{Mi^2 = I.4M0) as the sources detected 
by the ultimate DECIGO. We assume 
1 yr observation before the coalescence. We find that the signal-to-noise ratio is 
~ 2 X 10^(-ffo^S')~^) where Ds is the luminosity distance to the source. Fig|31 is 
the acceleration parameter X{z) in Eq. (|2-2j) as a function of z. The solid line is the 
case of ylCDM cosmology {h = 0.7, Qm = 1 — = 0.3). The filled circles and the 
error bars are the mean values of X{z) and the estimation errors {AX/Hq) /{ANY/"^ 
in Az = 0.1 bins. Thus, the cosmic acceleration rate at z = — 2 will be directly 
determined by 1 yr observation of the ultimate DECIGO. We note that the Hubble 
parameter H{z) is also measured by using Ea. (|2-2() : H[z) = (1 + z)[Hq — 2X[z)\ 
for the given Hq by another method. 

FigE] is similar to FigEl but the relative estimation errors {AX/ X) /{AN)^/'^ 
are shown for various observational period 1,3 and 10 yr. The errors are the smallest 
at 2; ~ 0.5 because X{z) is the largest there (see Fig|31). The errors are ~ 0.1 (10~^) 
for 1 (10) yr observation. In the matched filtering analysis, the phase modulation 
6\^acc in Ea. (|2-4() is roughly measured to an accuracy ^ {S/N)~^. Hence the relative 
error {AX/X)/{AN)^^'^ is estimated as 



AX 1 



X ./AN 



0.02 f ^/^° 



.05 



/ 



O.lHz 



At\ 
lyvj 



S/NY^ fANY^^^ 



10^ 



V 10^ 



(3-1) 

The above estimation is roughly consistent with the results in Fig0] The errors are 
proportional to At~'^ as shown in Ea. ()3-1() . For the high redshift z > 1 the error is 
large, since both S/N and X{z) become small. 

FigEl is the same as Fig01 but as a function of the binary mass with the fixed 
redshift of 2: = 1. The left panel is for the equal mass binary, and the right panel is for 
the unequal mass binary with the fixed companion masses Ml = 10^,10'' and IO^Mq. 
We assume that the event rate AN for the various mass binaries (10^"*^ — W^Mq) is 
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Fig. 3. The acceleration parameter X{z) (in 
Ea. (l2-2l l') as a function of redshift. The 
solid line is the case of ylCDM cosmology 
{h = 0.7, J7m = 1 ~ Qa = 0.3). The 
filled circles and the error bars are the 
mean values of X[z) and the estimation 
errors {AX/ Hq) / {ANf^ for NS/NS bina- 
ries (Mi,2 = IAMq) in Az = 0.1 bins. We 
evaluate for 1 yr observation before coales- 
cence by the ultimate DECIGO. 



Fig. 4. Similar to FigEl but the relative esti- 
mation errors {AX / X) / {ANY^'^ are shown 
for various observational period 1,3 and 10 

yr- 



the same as for the NS/NS binary for simpUcity. In the left panel, the errors are the 
smallest at a few M© mass, because the frequency at 1 — 10 yrs before the coalescence 
is near 0.1 Hz for these mass range and DECIGO/BBO is the most sensitive there. 
For a larger mass binary Mi_2 ^ lO^M©, though S/N is larger than that for NS/NS 
binary (1.4 -|- IAMq), the frequency is shifted to the lower one (< 0.1 Hz) so that 
the errors become large because DECIGO/BBO is less sensitive there. In the right 
panel, the observational period is fixed with Z\t = 3 yr but the result simply scales 
as (zAt/3yr)~^ from Ea. H3-l|) . For Mi = 10^ and 10*^ Mq, the errors are small for 
the smaller companion masses, and the reason is the same as in the left panel. For 
Ml = 10^ M0 with the smaller companion mass (M2 < lOM©), the S/N is small and 
hence the error is large. 

Finally, we discuss the estimation errors of the parameters of the dark energy 
by the ultimate DECIGO. The dark energy is usually characterized by its equation 
of state w(^z) = Pw{z)/ pw{z) where Pwiz) is the pressure and Pw{z) is the energy 
density. In the cosmology with the dark energy, the Hubble parameter at z is written 
as 

f l-\-w(z')M^^'^ 

H{z) = Ho + zf + (1 - - n^) (i + zf + n^e^ph dz' ' \ 

' (3-2) 

where is the density parameter of the dark energy. Note here that no assumption 
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Mi2(M®) M2(M®) 

Fie. 5. Same as FigEl but as a function of the binary mass with the fixed redshift z = 1. The left 
panel is shown for the equal mass binary, and the right panel is shown for the unequal mass 
binary (Mi — 10"^, 10^ and IO^A/q). In right panel, the observational period is fixed with 3 yr. 

on the curvature of the universe is made. We use the conventional first order expan- 
sion of w{z) as w{z) = wq+wiz where wq and wi are constants. Then the exponential 
in Eg . (13-21) is rewritten as -)-2)3(i+w)o--u)i)g3«)i2_ calculate the estimation errors 
of four parameters (i7m, ^w,W{),wi) instead of the acceleration parameter X{z) by 
inserting Ea. H3-2p into Ea. (|2-2() . We use the signal from all the NS/NS binaries at 
z = - 2 in FigEI 



obs. period AQm/Om AQw/^2w Awo Awi 



(a) 


w{z) — 


TOO 










1 yr 


0.030 


0.038 


0.031 






3 yr 


4.2 X 10"^ 


5.3 X 10"^ 


4.3 X 10"^ 






10 yr 


6.3 X 10"* 


7.6 X 10"'' 


5.9 X 10-"* 




(b) 


w{z) — 


Too + TOiz 










1 yr 


0.21 


0.62 


0.37 


0.81 




3 yr 


0.029 


0.086 


0.051 


0.11 




10 yr 


4.2 X 10"^ 


0.012 


7.3 X 10"^ 


0.016 



Table I. The estimation errors of cosmological parameters [fim, fi^^wo^wi) for various observa- 
tional period (1,3, 10 yr). The upper part is for (a) a constant to (i.e. w = too), and the lower 
part is for (b) a variable to (i.e. to = too + toiz). We use the sources of all the NS/NS binaries 
at z = - 2. 

In Tabled the results are shown for various observational periods i.e., 1,3 and 
10 yrs. The upper part is for (a) a constant w (i.e. w = t^o), and the lower 
part is for (b) a variable w (i.e. w = vjq + wiz). The results are shown for the 
cosmological parameters {Qm,^w,wq,wi) = (0.3,0.7,-1,0). From the upper part 
of Tabled these parameters {Qm, ^w,wq) are determined within a few percent by 1 
yr observation, which is comparable to the expected accuracy from SNAP satellite.^) 
We found that for the sources of small redshift {z < 0.5) the parameters are not well 
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obs. period 




Awo 






(a) w{z) = 


wo 








1 yr 


1.4 X 10"=* 


5.8 X 10" 


-3 




3 yr 


1.9 X 10"" 


8.0 X 10-^ 


-4 




10 yr 


2.9 X lO"'-' 


1.1 X 10" 


-4 




(b) w{z) = 


Wo + WlZ 








1 yr 


6.6 X 10"^ 


7.3 X 10" 


-3 


0.058 


3yr 


9.4 X 10"* 


9.6 X 10" 


-4 


8.0 X 10"^ 


10 yr 


1.4 X 10-* 


1.2 X 10' 


-4 


1.1 X 10"^ 



Table II. Same as TableEl but for flat universe model {Om + f^w = 1)- 

determined because these parameters are degenerate. But including the sources of 
high redshift (z > 1), the degeneracy is resolved and these cosmological parameters 
are well determined. 

Table ^ is the same as Table m but assuming the flat universe model (I7m + 
f2y^ = 1) suggested from the cosmic microwave background (CMB) observation by 
WMAP.io) The results are about 10 times better than that in Table 

We comment on the results for practical DECIGO or BBO. The errors by using 
these detectors are 1000 times larger than the results in the previous figures (FigsEl 
toEI) and tables (Table HI and |Tl|). But since the errors are proportional to At ^ {At is 
the observational period) from Ea. p-lf) . the cosmic acceleration could be measured 
with large observational period At > 10 yr. In fact, the parameters {f2m,WQ) will be 
determined within (3%, 10%) by 10 yr observation of practical DECIGO/BBO from 
Table llll(a). 

§4. Conclusion 

We have discussed the possibility of the direct measurement of the cosmic ac- 
celeration by using the decihertz gravitational wave detector. We take the binary 
of mass 0.1 - 10^ Mq at redshift z = - 2 as the source measured by the ultimate 
DECIGO. Following the previous paper, we newly include the effects of (1) the 
number of the sources and (2) the various masses and redshifts of the sources. We 
found that the expansion rate of the Universe (i.e. Hubble parameter) at 2; = — 2 
is directly determined within ~ 10% (0.1%) for 1 (10) yr observation by the ulti- 
mate DECIGO. The equation of state of the dark energy w is determined within 
~ 3% (0.06%) for 1 (10) yr observation. Even for practical DECIGO or BBO, w 
will be determined within ~ 10% by 10 yr observation assuming the flat universe 
model. Hence the DECIGO or BBO will determine the nature of the dark energy 
with high accuracy independently from the other cosmological tests such as SNAP, 
CMB observations. 
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